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Abstract 


■V 

Ws  adopt  tha  positioa  chat  spicttlea,  aacroapicalea  and  surges  ara 
■aaifastations  of  the  saaa  phanoaanon  occurring  on  different  scales.  We 
therefore  aearch  for  a  aechanisn  which  can  be  succesafnllj  applied  to  explain 
the  phenonenon  on  all  three  scales. 

We  first  consider  the  possibility  that  the  aechanisa  is  tha  same  as  that 
which  operates  in  producing  the  solar  wind*  except  that  Che  divergence  of  the 
aagnetic  ducts  is  auch  aore  rapid.  We  find  that  the  aechanisa  fails  to 
explain  spicules,  aacrospicules  or  surges.  For  instance,  if  it  produces 
speeds  typical  of  spicules,  the  aaxiaua  height  is  auch  too  aaall;  if  it 
reproduces  Che  height,  the  required  velocities  are  auch  too  high. 

We  also  consider  a  variant  of  this  aechanisa  proposed  by  Dcbida  in  which 
the  gas  pressure  is  suppleaented  by  the  aagnetic  pressure  of  a  gas  coaposed  of 
plasaoids.  This  aechanisa  also  fails  for  siailar  reasons. 


I.  I«t«o4— tioa 

Vhctt  th«  solar  liab  is  obsarvtd  with  s  narrow-band  filtar  eantarad  on  one 
of  Che  strong  chromospheric  emission  lines  (Ho  or  Ca  Z1  K),  the  chromosphere 
is  observed  to  consist  of  many  rapidly  changing  hairlike  features  called 
"spicules"  (Beckers  1972).  Each  spicule  appears  to  consist  of  a  tube  of 
diameter  about  1,000  km  extending  to  a  height  in  Che  range  6,000-10,000  km, 
with  a  velocity  in  the  range  20-30  km  s'^.  The  temperature  is  of  order 
10,000-20,000  K,  and  the  density  3.10^^  -  2^0^^em~^*  Bach  spicule  appears  to 
be  ejected  from  the  low  chromosphere,  and  it  may  fall  back  along  the  same  path 
and/or  fade  along  its  full  length.  The  total  lifetime  is  of  order  5-10  m. 
Spicules  occur  continually  on  the  surface  of  the  sun;  Beckers  (1972) 
estimates,  that,  at  any  time,  there  about  10^  spicules  on  the  sun's  surface. 
It  is  now  helieved,  on  the  basis  of  observations  of  spicules  seen  against  the 
solar  disk  (Simon  and  Leighton  1964)  that  spicules  tend  to  occur  at  the 
boundary  of  the  supergranulation  network.  It  is  also  worth  noting  that 
Sx  bright  mettles  tend  to  occur  at  the  roots  of  spicules  (Beckers  1972). 

In  a  projocted  seties  of  articles,  of  which  this  is  the  first,  we  plan  to 
examine  several  possible  mechanisms  for  spicules  and  also  for  surges  (Svestka 
1976;  Tandberg-Banssen  1974).  When  a  flare  starts,  it  is  quite  often 
accompanied  by  gas  ejections  which  manifest  themselves  as  moving  prominences 
on  the  limb  or  Doppler-shifted  dark  filaments  on  the  disk.  One  type  of  such 
ejection  is  a  "surge",  which  typically  has  the  form  of  a  straight  or  slightly 
curved  spike  which  grows  upward  from  the  chromosphere  with  velocities  in  the 
range  90-200  km  s'~^.  It  reaches  a  maximum  height  in  the  range  20,000-200,000 
km,  and  then  falls  back  to  the  chromosphere,  apparently  along  the  original 
trajectory.  The  lifetimes  of  these  spikes  are  in  the  range  10-30  m.  It 
appears  that  many  surges  are  preceded  by  a  diffuse  expansion  of  a  part  of  the 
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fl«r«.  toy  (1973)  has  idontificd  thoso  osponsioBs  with  "■ottotBchet"  or 
‘^lloraaa  boabo". 

Sorgo*  oottBlly  occur  at  aoall,  changing  aatellite  annspoti  which 
roproaant  ialanda  of  polarity  ravoraal  (tust  1968)  situated  close  to  the  edges 
of  sunspot  penunbrae  (Giowanelli  and  McCabe  1958).  When  satellite  sunspots 
disappear  froa  the  wicinity  of  a  spot*  surge  activity  ceases.  It  appears, 
therefore,  that  a  surge  is  due  to  the  aagnetic  field  reversal  of  a  satellite 
sunapot  and  that  the  process  which  gives  rise  to  a  surge  also  tends  to 
eliainate  the  field  reversal. 

Apart  froa  the  differences  in  scales  (length,  tiae  and  velocity),  there 
is  a  strong  siailarity  between  a  spicule  and  a  surge  (Tandberg-Hanssen  1974). 
If  a  surge  were  sealed  down  in  sise  and  in  velocity  to  valuea  typical  of 
spicules,  it  is  probable  that  it  would  appear  to  be  just  another  spicule. 

In  support  of  the  proposition  that  spicules  and  surges  are  due  to  basi¬ 
cally  the  saae  aeehanisa  occurring  on  different  scales,  we  also  point  out  that 
there  exist  'Siaero-spieules'*  which,  according  to  Beckers  (1977),  reseable  "a 
snail  surge"  or  a  "giant  spicule".  These  appear  in  the  polar  regions  of  the 
sun,  rising  to  heights  of  up  to  35,000  kn  with  velocities  of  up  to  150  kn  s*^, 
lifetimes  of  up  to  45  a,  and  diaaeters  of  up  to  10,000  ka.  Hence  the 
characteristic  diaensions  of  aacro-spicules  bridge  the  gap  between  those  of 
spicules  and  surges. 

For  reasons  indicated  briefly  above,  we  take  the  position  that  it  is 
reasonable  to  search  for  one  aeehanisa  responsible  for  spicules,  aacro- 
spicules  and  surges.  The  appropriate  values  of  physical  parameters  will  vary 
from  one  category  to  another,  but  it  is  reasonable  to  assume  that  the  saae 
basic  process  occurs  in  all  eases.  In  sesrehing  for  possible  aechanisas  of 
spicules  and  surges,  we  attach  priae  significance  to  the  fact  that  the  height 
of  the  spicule  or  surge  is  auch  greater  than  either  the  barometric  scale 
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height  or  the  hallistic  scale  height  for  the  relevant  velocity.  For 
teaperatores  usually  quoted,  the  haroaetric  scale  height  is  in  the  range  300- 
600  ka.  The  ballistic  scale  height  v^/2g  is  only  1,300  ka  for  a  typical 
spicule  velocity  of  25  ka  a~^,  and  is  only  20,000  ka  for  a  typical  surge 
velocity  of  100  ka  e'^.  It  is  clear,  therefore,  that  spicules  or  surges 
cannot  represent  siaple  ballistic  action  of  chroaospheric  aaterial  ejected 
froa  the  chroaosphere  by  one  aeans  or  enothsr.  This  aaterial  aust  be  snbject 
to  an  upward  force  which  counteracts  the  effect  of  gravity. 

It  is  also  significant  that  this  force  aust  be  coaparable  to  the 
gravitational  force.  Gravity  acting  alone  can  reduce  en  initial  flow  velocity 
of  20  ka  s~^  to  aero  in  a  distance  of  only  700  ka.  Siailarly,  if  the  upward 
force  were  (for  instance)  twice  that  of  gravity,  an  initial  flow  velocity  of 
20  ka  aec~^  would  be  increased  to  40  ka  aec”^  in  only  2,000  ka.  Hence  the 
spicule  aechanisa  auet  involve  a  force  which  is  coaparable  with  the  force  of 
gravity  and  reaains  fairly  constant  during  the  duration  of  a  spicule.  Siailar 
reaarks  apply  also  to  surges. 

In  this  article,  we  esaaine  briefly  two  aethods  of  providing  the  required 
auxiliary  force.  In  Section  II,  we  consider  the  possibility  that  a  spicule  is 
driven  by  the  pressure  of  the  spicule  gas  itself.  In  order  to  attain 
supersonic  speeds,  the  spicule  aust  then  be  produced  by  a  aechanisa  analogous 
to  that  responsible  for  the  solar  wind  (Parker  1963),  which  can  be  provided  by 
an  expanding  aagnetic  flux  tube.  In  order  to  siaplify  the  analysis,  we 
consider  en  idealised  aodel  in  which  the  flow  has  achieved  a  steady  state.  We 
find  that  this  aodel  cannot  produce  flow  heights  auch  greater  than  one  or  two 
ballistic  scale  heights. 

In  Section  III,  we  consider  a  proposal  of  Uehida  (1969),  which  aay  be 
regarded  as  a  aodification  of  the  above  aodel.  According  to  Dehida,  the 


Mttrial  floviat  upwards  io  a  spieala  wap  ba  eowposad  of  aaay  plaawa  bubblas 
or  "plasaoida".  la  this  caao,  the  •ataatle  praasoro  of  tba  plasaoid 
tapploaaacs  ebo  (as  prossuro.  Vo  follow  Uehida  la  coasidorlag  a  stoady-state 
flow:  this  iawolras  oaly  a  slight  oxtaasioa  of  tbo  aaalysis  of  Soetioa  II. 
Oaeo  agaia,  wo  fiad  that  this  aodol  doos  aot  prodoeo  flow  haights  wueb  graatcr 
tbaa  oaa  or  two  ballistic  seals  haights. 

Tba  failara  io  aaeb  of  thsss  modals  to  produea  a  largo  ratio  of  flow 
boigbt  to  ballistic  height  may  bo  traced  to  the  diffieolty  of  prodociog  ao 
almost  eoastaat  force  io  ao  eapaadiag  regioa.  la  the  first  model,  the 
espaasioo  which  is  accessary  for  the  aeeeleratioa  also  causes  ceoliog  of  the 
spicule  material.  The  aeeeleratioa  is  doe  to  pressure,  but  the  pressure  is 
reduced  by  eapaosioo  aod  the  force  becomes  oegligible  as  sooo  as  the  gas 
cools.  Io  the  seeoad  model,  the  magaetic  force  depeads  seositively  oo  the 
area  of  the  tube.  If  large  eapaosioo  occurs  withia  the  flow  regioa,  the 
magaetic  force  will  aot  remaio  comparable  to  the  grawitatioaal  force,  la  this 
ease  either  the  spicule  velocity  becomes  too  large  or  the  spicule  height 
becomes  too  small.  If  the  scale-height  for  eapaosioo  of  the  tube  is 
iacreased,  the  effects  of  eapaosioo  are  reduced,  but  this  reductioa  cao  be 
offset  by  aa  iacrease  io  the  assumed  magaetic  field  streogth  ia  the  plasmoids. 
la  this  way,  the  parameters  io  the  Uchida  model  cao  iadeed  be  adjusted  so  that 
the  magaetic  force  oearly  caocels  gravity  aloog  the  length  of  the  spicule. 
However,  this  solution  is  very  sensitive  to  fluctuations  in  the  flow  velocity, 
la  the  absence  of  a  neebanisn  leading  to  this  delicate  adjustment  of 
parameters,  the  Uchida  model  does  aot  seem  to  explain  spicule  behavior. 

II.  Steady-State  Hydrodvamaic  flow 

He  consider  steady  flow  along  a  thin  flux  tube  which  begins  and  returns 
to  the  chromosphere  which  thus  provides  both  a  source  and  a  sink  for  the  gas 
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is  th«  spicttl*.  However,  einee  Che  beaie  requireweat  of  «  spicule  model  is 
that  it  explaia  the  upward  motion  of  chromospheric  gas.  we  shall  restrict  our 
atteatioa  to  the  initial  upward  motion  and  ignore  the  subsequent  downward 
flow. 

For  simplicity,  we  consider  steady  upward  flow  along  a  thin  vertical  flux 
tube  of  which  the  area  A  is  a  given  function  A(s)  of  a.  where  s  is  measured  in 
the  upward  direction.  Then  the  equation  of  continuity  is 

nvA  -  J  .  (2.1) 

where  n  is  the  particle  number  density,  v  the  velocity,  and  J  (the  particle 
mass  flux  along  the  tube)  is  independent  of  s.  All  quantities  are  measured  in 
e.g.s.  units. 

The  equation  of  notion  is 


nnv 


dv 

ds 


(2.2) 


where  m  is  the  mean  particle  mass,  g  the  (downward)  gravitational  acceleration 
and  p  the  gas  pressure. 

The  energy  equation  may  be  written  as 

3  dT  p  d 

2  dJ  “  ■  A  d7  ^  »  <*-5) 


where  k  is  Boltsmann's  constant.  T  is  the  temperature,  and  Q  is  the  net  rate 
of  energy  input  per  unit  volume,  taking  account  of  energy  loss  by  rsdiation. 
For  the  low  temperatures  and  other  conditions  characteristic  of  spicules  and 
surges,  the  role  of  thermal  conduction  in  the  energy  equation  may  be  ignored. 
To  the  above  equations  we  must  add  the  relation 

p  -  nkl  .  (2.4) 
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For  tho  toaporotorts  and  oChor  eonditioiia  of  iottrtst,  ionisatioo  will  play 
only  a  aaall  zolo  and  will  be  noslaetad.  ffa  tharefora  conaidar  tha  aaaa  a  to 
ba  tha  aaaa  of  a  hydrogen  atoa.  On  tha  other  hand,  va  conaidar  that  there  ia 
aofficiant  ioniaation  to  inaure  that  flow  ia  along  aagnetie  field  linea. 

We  nay  naa  aquation  (2.4)  to  alininata  p,  aquation  (2.1)  to  eliminate  n, 
and  then* rearrange  aquationa  (2.2)  and  (2.3)  to  obtain  espreaaiona  for  dv/dc 
and  dl/da.  Tha  reaulting  aquationa  are  aa  followa: 


/  c3\  dw  .  C*  dA 

2  AQ 

V  T  /  da  *  '  A  da 

3  Jn 

1  "  -  •  -  i  “  ♦! 

f5 

1 

\  vV  T  da  ^ 

^3 

*  JaC^ 

where  C  ia  the  sound  speed  defined  by 


(2.5) 


(2.6) 


(2.7) 


We  see  that  each  of  these  equations  has  a  "critical  point"  defined  by  the 
saae  conditions: 


(2.8) 


and 


Ao  Cj  3 


(2.9) 


For  given  physical  conditions,  such  as  flux-tube  geonetry,  heat  input  and 
radiation  losses,  one  nay  view  equations  (2.8)  and  (2.9)  as  representing  a 


regulating  nechanisa  for  thi 


s  flux  J. 


For  purpoMS  of  eoaputiog  lolutiono  of  tho  above  equetiooSf  it  le 
cosveaient  to  expxeee  thea  la  diaensionleoe  fora.  Ve  therefore  introduee  the 
folloviag  diaeneioalees  verieblee. 


C  •  */®o 
•  -  A/A^ 

€  -  v/C^ 

T  -  T/I^ 

U  ■  Q/Qo 


end  the  following  diaeneionlees  paraaetere. 


r 

A 


cB  C“^ 
•“o'*© 

2  Vo'S. 
5  .JC* 


In  thcae  expreaaiona,  ia  defined  by 


(2.10) 


(2.11) 


(2.12) 


The  quantity  F  ia  the  ratio  of  the  scale  length  characteristic  of  the 
expansion  of  the  flux  tube  at  the  critical  point  to  the  ballistic  scale  height 
at  the  critical  point.  The  quantity  A  aay  be  viewed  approxiaately  as  the 
ratio  of  the  rate  of  energy  input  in  the  neighborhood  of  the  criteal  point  to 
the  kinetic  energy  flux  at  the  critical  point. 

To  these  expressions,  it  is  convenient  to  introduce  a  syabol  denoting  the 
Mach  nuaber  of  the  flow: 


N  « 


V 

C 


(2.13) 
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lo  terms  of  dimeas  ion  lest  vsrisbles,  equations  (2.5)  and  (2.6)  now  become 


(1  -  M"*)e 


d( 


1  da 

-  r  ♦  ;  -  -  A.U 


(2.14) 


and 


I  (1  .  H->)  £  -r^  -  i  il,,  I  (|  -  .  <2.») 


We  see  from  equations  (2.8)  and  (2.9)  that  the  conditions  at  the  critical 
point  may  now  be  expressed  as 


-  1  (2.16) 

and 

r  ♦  A  -  1  .  (2.17) 

It  follows  from  (2.17)  that  T  <  1,  i.e.,  that  the  expansion  scale  height 
less  than  cjg  which  is  twice  the  ballistic  scale  height 
corresponding  to  the  speed  of  sound  at  the  critical  point.  It  also  follows 
(2.17)  that  A  <1.  (We  are  assuming  that  there  is  a  net  energy  input 
at  the  throat  so  that  A  >  o.) 

We  have  computed  solutions  of  equations  (2.14)  and  (2.15)  for  a  range  of 

assumed  forms  for  the  area  function  A(s)  and  heating  function  Q(x).  A  typical 

solution  is  shown  in  Figure  1.  For  this  particular  case,  the  area  A(s)  is 
assumed  to  grow  exponentially  from  the  throat  on,  as  follows:- 

10  for  s  <  0  • 

l/(  500  km)  for  0  ^  a  <  2000  km  ,  (2.18) 

1/(5000  km)  for  s  ^  2000  km  . 
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For  the  heeCing  functioa,  «e  adopt  the  sauaeian  form 


Q(z) 


Qo  «*P 


I 


(2.19) 


which  peaka  at  the  value  Xp  and  haa  a  width  determined  by 

The  velocity  v(z)  and  temperature  T(z)  reaulting  from  thia  model  are 
ahown  in  Figure  1.  It  ia  to  be  noted  that  the  temperature  decreaaea  very 
rapidly,  due  to  adiabatic  expanaion.  The  velocity  increaaea  for  temperaturea 
greater  than  10,000  K  but  decreaaea  rapidly  when  the  temperature  dropa  below 
10,000  K.  It  ia  clear  that  thia  particular  model  doea  not  meet  the  require¬ 
ment  of  a  apicule  model,  that  the  velocity  remaina  conatant  over  aeveral 
balliatic  acale  heighta. 

The  above  reault  ia  not  peculiar  to  the  caae  apecified  by  equationa 
(2.18)  and  (2.19);  it  aeema  to  be  quite  general.  Thia  property  of  the 
hydrodynamic  model  may  be  underatood  aa  followa.  In  order  for  the  flow  to 
continue  to  accelerate  after  it  becomea  auperaonic,  the  equation  of  motion 
(2.14)  muat  be  dominated  either  by  the  expanaion  term  (1/a)  (da/d^)  becoming 
large  and  poaitive,  or  by  the  heating  term  AaM  becoming  large  and  negative 
(repreaenting  a  large  heat  loaa  by  radiation).  In  either  caae,  we  aee  from 
equation  (2.15)  that  the  effect  ia  to  lead  to  a  negative  value  of  dr/d^  ,  ao 
that  the  temperature  decreaaea  with  height.  However,  if  the  temperature 
decreaaea  with  height,  the  expanaion  term  in  equation  (2.14)  will  aoon  become 
negligible.  The  flow  then  becomea  eaaentially  balliatic.  Balliatic  flow  and 
a  decreaae  in  temperature  with  height  are  both  inconaiatent  with  apicule 
obaervationa  (Beckera  1972). 

Thia  model  alao  faila  to  reproduce  the  propertiea  of  aurgea.  In  order  to 
attain  heighta  characteriatic  of  the  largeat  aurgea  (~  10^^  cm),  we  require 
that  Cg/2g  attain  thia  value.  Thia  then  requirea  that  >  10^*^,  which  in 
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turn  requires  e  teapereture  T  :  10^*^.  This  teapersture  is  aucb  higher  Chen 
the  ceapereture  of  surge  aster iel,  which  is  of  order  10^  K. 


III.  Modified  Pehide  Model 

Dehids  (1969)  has  considered  a  aodel  which  is  a  variant  of  that  con¬ 
sidered  in  Section  II.  He  considers  that  plasaoids  are  produced  during  the 
reconnection  process,  each  plasaoid  containing  gas  and  a  closed  aagnetic 
field.  The  gas  pressure  is  Chen  suppleaented  by  aagnetic  pressure. 

In  this  section,  we  follow  Uchida  in  this  general  concept.  However, 
Uchida  assuaes  that  the  aagnetic  field  strength  of  the  plasaoids  is  a  known 
function  of  height.  Ve  prefer  to  estiaate  the  aagnetic  field  strength  of  the 
plasaoids  by  an  "equation  of  state".  Ve  asauae  that  the  aagnetic  field  is 
"frozen"  into  the  plasaa.  Then  the  aanner  in  which  the  aagnetic  field 
changes,  in  response  to  a  change  in  the  gas,  depends  on  what  assuaptions  one 
aakes  concerning  the  relative  directions  of  expansion  and  of  the  aagnetic 
field.  Tor  instance,  if  the  expansion  is  along  the  aagnetic  field,  then  B  is 
independent  of  the  gas  density  n.  If,  on  the  other  hand,  Che  expansion  is 
transverse  to  the  aagnetic  field,  then  B  is  proportional  to  n. 

He  assume  that  the  magnetic  field  is  statistically  isotropic,  from  which 
it  follows  that  B  «  n^^^.  Since  the  magnetic  pressure  is  proportional  Co  B^, 
it  follows  that 


He  nay  therefore  write 


(3.2) 
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where  end  will  be  interpreted  ••  the  veluee  at  the  "throat".  Od  uaing 
equations  (2.1)  and  (2.10),  (3.2)  hecoaes 

"  Pa.o  5  •  <3-3) 

The  analysis  of  this  section  follows  closely  that  of  the  previous 
section.  The  equation  of  continuity  (2.1)  and  the  energy  equation  (2.3)  are 
unchanged.  However,  the  equation  of  notion  (2.2)  is  now  replaced  by 


dv 

ds 


dp  dp^ 
ds  ds 


(3.4) 


Another  effect  of  the  nagnetic  pressure  is  to  change  the  wave-propagation 
speed,  so  the  speed  of  sound  C  is  now  replaced  by  K,  where 


,2  5  KT  4  H 

*  "3r*3« 


(3.5) 


The  "Mach  nuaber"  will  now  be  defined  as 


"-i 


(3.6) 


In  equations  (2.10)  and  (2.11),  is  now  replaced  by  K^,  so  that  the 
diaensiottless  variables  and  paraaeters  becone 


C  -  */Ho 
a  -  A/Ap 
C  -  v/K^ 

T  -  T/T^ 

P  -  Q/Qo 


(3.7) 
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and 


It  is  convenient  nleo  to  introduce  the  symbol 

a  -  , 


(3.8) 


(3.9) 


which  may  be  viewed  as  a  normalised  "effective  temperature". 

With  this  terminology*  we  now  obtain*  in  place  of  (2.14)  and  (2.15)*  the 
equations 


(l-ir*)e^--r  ♦  o-~-Aay 


(3.10) 


and 


Ve  see  that  this  pair  of  equations  also  allows  a  transition  from  subsonic  to 
supersonic  flow  at  a  critical  point  at  which  equations  (2.16)  and  (2.17)  most 
be  satisfied. 

Two  numerical  solutions  of  these  equations  are  shown  in  Figures  2  and  3, 
the  former  representing  a  case  with  moderate  magnetic  field  strength*  and  the 
latter  a  ease  with  strong  magnetic  field  strength.  The  magnetic  stress  leads 
to  higher  peak  velocities  and  greater  heights  of  the  flow  patterns.  As 
before*  rapid  cooling  occurs  in  the  expansion  region*  and  the  motion  is 
essentially  ballistic  beyond  the  height  at  which  the  velocity  has  its  maximum 
value.  In  the  ease  studied  in  Figure  2*  the  maximum  velocity  (50  km  a"^)  is 
larger  than  ia  typical  of  spicules*  yet  the  model  does  not  extend  to  heights 
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typical  of  spieulaa.  la  the  eaae  atudlad  ia  Figura  3,  the  flow  catcads  to 
hoighta  typical  of  apieulaa,  but  tha  aaaiaua  oalocity  (75  ka  b~^)  ia  much 
larger  thaa  valoaa  typical  of  apieulaa.  Froa  thaaa  and  other  axaaplaa,  it  ia 
clear  that  the  Uchida  aodel  doea  aot  aeet  the  requireaeata  of  a  valid  aodel  of 
apieulaa . 

Ve  have  alao  iaveatigated  the  poaaibility  that  the  Dehida  aodel  aay 
repreaeat  the  aeehaaiaa  of  aurgea.  However,  the  aodel  faila  for  the  aaae 
reaaoaa  that  it  faila  to  eaplaia  apieulaa.  If  the  apeed  ia  correct,  the 
height  ia  too  aaall;  if  the  height  ia  correct,  the  apeed  ia  too  high. 

If.  Biacaaaioa 

Va  have  ahowa  that  two  of  the  eoaeeptually  aiapleat  aodela  of  apieulaa 
aad  aurgea  fail  to  reproduce  their  propertiea.  loth  phettoaeua  require  that 
the  gaa  aove  under  the  influence  of  a  force  which  alnoat  balancea  the  force  of 
gravity  over  a  height  which  ia  nuch  larger  thaa  the  baroaetric  aeale  height 
aad  aueh  larger  than  the  balliatie  acale  height.  Heither  of  the  aodela 
providea  auch  a  force. 

Even  if  the  above  iM»dela  had  been  able  to  reproduce  the  velocity  enrvea 
of  apiculea  and  aurgea,  they  would  have  encountered  another  difficulty.  In 
order  to  accelerate  gaa  in  a  duct  to  auperaonic  apeeda,  it  ia  neceaaary  that 
the  duct  expand.  The  tranaverae  dinenaion  of  the  gaa  would  therefore  expand 
with  height.  However,  obaervational  data  indicate  that  the  tranaverae 
dinenaion  of  apiculea,  if  anything,  decreaaea  with  height  (Beckera  1972). 

There  are  aeveral  varianta  of  the  above  nodela  which  can  be  conaidered. 
It  ia,  for  inatance,  poaaible  that  the  aaaunption  of  ateady  flow  aonehow 
rulea  out  an  eaaential  aapect  of  the  apieule/aurge  nechanian.  It  ia  poaaible 
that  the  cool  gaa  which  ia  viaible  in  a  apicule  or  aurge  ia  really  being 
expelled  by  the  preaaure  of  a  very  hot  gaa,  juat  aa  a  low-tenperature  bullet 
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is  sxpsllsd  by  hitb-trapcrstars  gas  ia  tbs  bsrrtl  of  s  rifls.  Boworor,  thoss 
▼sriaats  art  likaly  to  faeo  tba  saaa  objoetioas  as  aotod  abovs  eoaesraiag  tho 
tzaasTsrsa  diaoasioa. 

Aaotbar  possibility  is  that  tba  driviag  fores  for  a  spicule  or  surge  ie 
pi^ouided  by  aagaetic  teasioa  rather  tbaa  gas  pressure  or  aagaetie  pressure. 
This  suggestioa  is  embodied  ia  a  proposal  adraaead  some  time  ago  by  Fikel'aer 
(1969).  This  possibility  appears  the  most  proaisiag,  sad  will  be  discussed  ia 
our  ueat  article. 

This  work  was  supported  ia  part  by  Mkgk  Create  MGL  05-020-272  sad  I&6W-92 
aad  the  Office  of  Mawal  Beseareh  Coatraet  I00014-75-C-0673. 
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neon  c&mow 


Figure  1 


Figure 


Figure  3 


Hydrodynemic  model  with  the  critical  temperature  equal  to  20,000 
degreea.  The  ezpanaion  and  beatiog  funetiona  are  ahovn  in  Figure 
la).  The  velocity,  ahovn  aa  a  function  of  height  in  lb),  decaya  in  a 
diatance  comparable  to  the  balliatic  acale  height  of  the  maximum 
velocity.  The  temperature  alao  ia  ahovn  in  Figure  lb);  note  the 
rapid  cooling  that  occura  due  to  the  large  expanaion  near  the 
throat. 

Oehida  model  vith  moderate  magnetic  field  (V^>10  km  a**^,  C0-25  km  a~^ 
at  the  critical  point).  Comparing  thia  vith  Figure  1,  ve  find  a 
larger  peak  velocity  and  apicule  height  due  to  the  magnetic  force, 
▲a  in  the  hydrodynamic  model,  rapid  cooling  occura  in  the  expanaion 
region,  and  the  motion  ia  eaaentially  balliatic  paat  the  peak 
velocity. 

Uchida  model  vith  atrong  magnetic  field  (7^"20  km  a~^,  Cq*25  km  a~^). 
The  acceleration  occura  in  the  region  of  rapid  expanaion.  The 
apicule  height  ia  not  much  greater  than  one  balliatic  acale  height 
paat  the  peak  velocity.  While  the  height  ia  cbaracteriatic  of  long 
apicttlea,  the  peak  velocity  (7$  km  a'^)  ia  much  larger  than  obeerved 
apicule  velocitiea. 
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